00 

o 

O 
(N 



C/3 
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Multi- valley spin relaxation in n-type GaAs quantum wells with in-plane electric field is investi- 
gated at high temperature by means of kinetic spin Bloch equation approach. The spin relaxation 
time first increases and then decreases with electric field, especially when the temperature is rela- 
tively low. We show that L valleys play the role of a "drain" of the total spin polarization due to the 
large spin-orbit coupling in L valleys and the strong T-L inter- valley scattering, and thus can enhance 
spin relaxation of the total system effectively when the in-plane electric field is high. Under electric 
field, spin precession resulting from the electric-field-induced magnetic field is observed. Meanwhile, 
due to the strong T-L inter-valley scattering as well as the strong inhomogeneous broadening in L 
valleys, electron spins in L valleys possess almost the same damping rate and precession frequency 
as those in V valley. This feature still holds when a finite static magnetic field is applied in Voigt 
configuration, despite that the ^-factor of L valleys is much larger than that of F valley. Moreover, 
it is shown that the property of spin precession of the whole system is dominated by electrons in 
r valley. Temperature, magnetic field, and impurity can affect spin relaxation in low electric field 
regime. However, they are shown to have marginal influence in high electric field regime. 
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I. INTRODUCTION 

The ability to manipulate electron spin degree of free- 
dom in semiconductors has become an important issue 
for the realization of spintronic devices.— i^i^ Usual ways to 
manipulate spins arc proposed by directly applying mag- 
netic field or exploiting the spin-orbit interaction through 
a gate voltage^i^ and/or strain field i^'^'^'iS In addition, 
applying a drift electric field also provides an approach to 
manipulate spin, by which the dependence of spin-orbit 
coupling on electron momentum is exploited and an ef- 
fective magnetic field can be inducedjii However, these 
manipulations also affect spin relaxation and dephasing 
(R&D) time,^-i^ii^'i^ii^ which is an essential time scale 
for the design of spin-based devices. Moreover, as the 
goal of semiconductor spintronics is to combine "tradi- 
tional" semiconductor electronics with the utilization of 
the spin state and most current electronic devices are 
performed at high-field transport condition, the study of 
drift electric field on spin R&D is therefore essential. 

Earlier Monte-Carlo simulation has revealed that 
the drift electric field can enhance spin relaxation in 
GaAs quantum wire a^^i^^ and bulk material f^^ How- 
ever, these studies fail to treat all scattering, especially 
the Coulomb scattering, completely. It was pointed 
out that the Coulomb scattering can also cause spin 
R&D.—'— Investigations based on fully microscopic ki- 
netic spin Bloch equation (KSBE) approach^^'^^ have 
demonstrated that the Coulomb scattering is very im- 
portant to the spin R&D i^^i^^i^^'^° i^ This has also been 
verified experimentally.— i^^i^l Among these works, Ref. 
Ill] gives a complete understanding of the hot-electron ef- 
fect on spin R&D in n-type GaAs quantum wells (QWs), 



where all the scattering is explicitly included and calcu- 
lated self-consistently. It is shown that the spin R&D 
increase with electric field at high temperature, but de- 
crease with it at low temperature^^ Nevertheless, due 
to the lowest-subband used in the investigation and the 
"runaway effect" ^^ that study is limited to electric field 
lower than 1 kV/cmiii When the electric field is fur- 
ther raised, the multi-subband or even the multi-valley 
effect should be taken into account i^ In Ref. [20, the 
electric field has been raised up to 3 kV/cm and the 
multi-subband scattering is considered. It is shown that 
the spin R&D first decrease then increase with the ap- 
plied electric field at low temperature.— These effects^i^S 
have very recently been demonstrated by HoUeitner et 
al. experimentally>S§. However, when the electric field is 
further increased, the inter-valley scattering should be 
taken into account. 

The inter-valley scattering can lead to interesting phe- 
nomena, such as the charge Gunn effect in GaAsjSI which 
comes from the different mobilities in different valleys 
due to the different effective masses. Studies focusing on 
scattering and transport between different valleys with- 
out the spin degree of freedom have been performed long 
time ago.— i^i^SiSii^^ When spin is considered, the dif- 
ferent spin-orbit couplings and the different effective g- 
factors, together with the well known different momen- 
tum relaxations at different valleys may lead to complex 
spin R&D phenomena. Spin injection through a Schot- 
tky barrier into bulk GaAs, with T-L-X valley structure 
included, has been studied by Monte-Carlo simulation, 
revealing faster decay of current (electron) spin polariza- 
tion in the upper valleys.— Recently, a spin Gunn effect 
has been proposed, depicting a process of spontaneous 



spin-polarization amplification in semiconductors such as 
GaAs Ji However, as this effect is predicted to happen in 
the high-electric-field charge Gunn domain, and the con- 
dition for it depends strongly on spin relaxation time, 
it is necessary to investigate spin relaxation under high 
electric field. 

In this work, we perform a fully microscopic investiga- 
tion on spin relaxation with high in-plane electric field 
at high temperature in n-type GaAs QWs, with F and 
L valleys included. The spin relaxation comes from the 
D'yakonov-Perel' mechanisnii^ This paper is organized 
as follows. In Sec. II, we set up the model and construct 
the KSBEs. We present our main results in Sec. Ill, in- 
cluding the role of inter- valley scattering, the spin preces- 
sion properties, and the effects of temperature, magnetic 
field, and impurity on spin R&D. We summarize in Sec. 
IV. 



II. MODEL AND KSBES 

We start our investigation from an 7i-type GaAs [001] 
QW with well width a. The growth direction is as- 
sumed along z-axis. F valley is located at the cen- 
ter of the Brillouin zone and four L valleys are at 
K" . = (7r/ao)(l, ±1, ±1) that are rotationally symmetric 
around the [100] {x) axis, with ao denoting the lattice 
constant and i = 1,- • • ,4. Under spherically symmet- 
ric approximation, the effective electron masses of F and 
L valleys are 7TT-r=0.067mo and 771^^ =771^=0. 230mor^'^ 
where toq is the free-electron mass. Assuming parabolic 
and isotropic band structure for each valley, the electron 
energy at the bottom of each valley reads^ 
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where k^ is the x-y-plane projection of the three- 
dimensional relative momentum K — K^ (A=r, Li), with 
the electron energy at Kp set as the reference point. 
£rL is the effective T-L energy difference in the two- 
dimensional structure: epL — Eyl — %^{— —), 

where £'rL=0.28 eV is the T-L energy difference in 
bulk;^ Here an infinite-well-depth assumption is used. 
a is set as 7.5 nm throughout the paper, etl — 0.21 
eV is smaller than the energy difference of the lowest 
two subbands of the F valley, which is about 0.30 eV. 
It is also noted that the effective L-X energy difference 
in this two dimensional structure is as high as 0.17 eV 
when the spherical effective mass of X valleys is about 
O.Bttiqj^ Thus in this work, only the lowest subbands of 
F and L valleys are considered. 

The Hamiltonian of the system reads 



^ = E^"+E^Av 



AA' 



H 



XX' 



Hrx'+Hl7 + HrySxx' 



(3) 
(4) 



H", 



^ {(ffk, +eE-R)5, 



kxaa' 



[-gx^iBB + h)^ik^)] ■ o-^^,}cLa 



Ckx'y' 



(5) 



is the free electron Hamiltonian at A valley, where — e 
is the electron charge, R=(a;, y) is the position of the 
electron, E is the electric field (the system is assumed to 
be spacial homogeneous one so that E does not depend 
on the position), B is the magnetic field in Voigt config- 
uration and <T is Pauli matrices. In the calculation, E 
and B are taken along x-axis unless otherwise specified. 
Cj^ cr(^kAO-) is the creation (annihilation) operator for 
electron with relative wave vector k^ and spin a. The ef- 
fective (/-factors are gr = —0.44^^ and gLt = gh ~ 1.77J^ 
hA(k) represents the Dresselhaus spin-orbit coupling4£ 
In the coordinate system adopted in this work^^ 
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Here n^ is the unit vector along the longitudinal prin- 
ciple axis of Li valley. {k\z) {{kf^)) represents the av- 
erage of the operator —id/dz — K^^ [{—id/dz — K\^'^\ 
over the electron state of the lowest subband in A val- 
ley. With the infinite-depth assumption, it reads (fc>,z)=0 
[(fc^^)=(7r/a)^]. The spin-orbit couphng coefficients uti- 
lized in the calculation are 7 = 0.011 eV-nm'^~'™ and 
/? = 0.026 eV-nm^ respectively. The interaction Hamil- 
tonian is composed of the intra- (A' = A) and inter-valley 
(A' 7^ A) Coulomb scattering iJ^^,^, electron-phonon scat- 
tering -ff^^/'' and electron-impurity scattering H'^f . Ex- 
pressions of these Hamiltonian can be found in books^i^ 
and Refs. 28,46. 

The KSBEs constructed by the nonequilibrium Green's 
function method read^^ii^ 



P^x 



P^\\dri + P'^\\coh + 



7 . P^x\scat,XX' 
A' 



(8) 



Here pk^ represent the density matrices of electrons with 
relative momentum k^ in valley A, whose diagonal terms 
PkA,CTCT = /k^.o- (cr = ±1/2) represent the electron dis- 
tribution functions and the off-diagonal ones Pt 1 1 = 
p* 1 1 describe the inter-spin-band correlations for the 

•^A, 22 

spin coherence. PkAJdri = eE • Vk^PkA are the driving 
terms from the external electric field. pkA \coh are the co- 
herent terms describing the coherent spin precessions and 
PkA [scat, AA' stand for the intra- (A = A') and inter-valley 
(A ^ A') scattering terms. Expressions of these terms are 
given in Appendix A. 

The initial conditions at time t—Q are prepared from 
spin unpolarized equilibrium states (in the absence of 
electric field) at time t — —to: 

hrA-to) = {exp[(£kr - M)/fcsr] + l}'^ , (9) 

h,„A-to) - Pkr.i-i(-^o) = Pk, i-i(-io) - (10) 



with /i representing the chemical potential at tempera- 
ture T. We turn on the electric field at i = — tg and 
allow the system to evolve to the steady state before 
i = 0. Under the driving field and the inter- valley 
scattering, part of electrons in F valley are driven to 
L valleys. Then we turn on a circularly polarized laser 
pulse with width 5^ = 0.01 ps at i = to excite spin- 
up electron with Gaussian-like distribution to T valley: 
'^/kr,i/2=aexp[-(ekr - £f)V2'^?][1 - /kr, 1/2(0)]- Here 
a = iVFo/IEkr exp[-(ekr " epfmW^ " /kr,i/2(0)]}, 
ep is the Fermi energy of F valley and 5ir=h/5r- N is the 
total electron density after the pulse excitation and Pq is 
the spin polarization excited by the pulse. In this work, 
TV = 4 X 10" cm-2 and Pq = 5 %. 




FIG. 1: (Color online) A typical time evolution of the ini- 
tial steady-state preparation (density A''a and spin polariza- 
tion Pa) under the electric field P = 4 kV/cm and a spin- 
polarization generation by a laser pulse at i = 0. Ni — B = 0. 
Chain curves: F valley; Dashed curves: L valleys; Solid 
curves: the total. 



ential electric conductance of the total drift velocity is 
obtained. 

The hot-electron temperature Te can be obtained by 
fitting the Boltzmann tail of the calculated steady-state 
electron distribution of each valley. ^^ Our results show 
that each valley has its own hot-electron temperature(s). 
Moreover, for F valley, there are two temperatures. One 
is at the energy regime which overlaps with that of L 
valleys. In this regime, electrons at F and L valleys share 
the same hot-electron temperature (labeled as Tl) due 
to the strong inter-valley scattering. The other is Tp for 
electrons in the lower energy regime of F valley. That is 
the typical hot-electron temperature of F valley and is 
higher than T^. Tr and Tl as functions of electric field 
at different lattice temperatures are plotted in Fig. [3l It 
shows that electrons with smaller effective mass (such as 
those of F valley) and lower lattice temperature (thus 
relatively weak scattering) are easier to be accelerated 
and heated. 
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FIG. 2: (Color online) Steady-state drift velocity v\ against 
electric field P at T=100 K. iVi = P = 0. Chain curve: F 
valley; Dashed curve: L valleys; Solid curve: the total. 



III. NUMERICAL RESULTS 

We numerically solve the KSBEs following the scheme 
mainly laid out in Ref. I 111 , with extension to the inter- 
valley scattering given in Appendix B. In Fig. [T] we 
show a typical time evolution of the initial steady- 
state preparation under the influence of the electric 
field and a spin-polarization generation by a laser pulse. 
The electric field is applied from time t — —to — 
—5.5 ps, and the laser pulse is applied at t = ps. 
Pxit) ^ 2Ek,.[^/k...(<)]/Ek../k.,.W and Nxit) ^ 
X^k a fkx,ui't) are the spin polarization and the electron 
density at A valley, separately. The steady-state drift ve- 
locity of A valley va can be obtained from the steady 
value of VA(i) = 'Ek^Ahx.<yit)lT'^\/'m\]/ J2k^a hx,<T{t). 
The steady-state drift velocities v\ of each valley as well 
as the total drift velocity against the electric field E are 
plotted in Fig. [21 From the figure the negative differ- 



The spin relaxation time r can be obtained from the 
temporal evolution of spin polarization Pa. Due to the 
lowest conduction subbands used in each valley and the 
"run-away" effect, ^'^ our research is limited to electric 
field up to 7 kV/cm. 



A. Temporal evolution of spin polarization 

We first show the temporal evolutions of total spin po- 
larization Protaiit) = J2\P^i^) under different electric 
fields at P = 200 K with P = iV, = in Fig. H The 
corresponding spin relaxation times obtained from the 
calculations are also given in the figure caption. It is 
seen from the figure that in the presence of electric field, 
electron spins precess with the precession period show- 
ing P-dependence. This is because the center-of-mass 
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FIG. 3: (Color online) Hot-electron temperature Tr (curves 
with squares) and Tl (curves with circles) against electric 
field E at lattice temperature r=300 K (chain curves) and 
200 K (dashed curves). iVi = B = 0. 



effect is quite marginal compared to the first oneiii It 
has been demonstrated before that electrons in F valley 
is in the strong scattering limit jiii^idl When the elec- 
tric field is low, electrons are mainly distributed in the 
F valley. Thus the influence of hot-electron effect in the 
strong scattering limit causes r to increase with E. How- 
ever, when the electric field is high enough, apart from 
the fact that the increase of inhomogeneous broadening 
from the cubic term of the Dresselhaus term in F val- 
ley becomes important)^ L valleys also start to play an 
important role. The amount of electrons sitting in L 
valleys grows with the electric field and the hot-electron 
temperature increases with the field as well. Thus the 
inter-valley scattering becomes more important with the 
increase of electric field. The strong spin-orbit coefficient 
at L valleys in conjunction with the inter-valley scatter- 
ing lead to the decrease of the spin relaxation time with 
the field. 



B. Effect of inter- valley scattering on spin 
relaxation 
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FIG. 4: (Color online) Temporal evolution of spin polariza- 
tion Protai when E = (solid curve), 2 (dashed curve), and 
4 kV/cm (chain curve). The corresponding spin relaxation 
times are fitted to be 10.5, 14.4, and 4.7 ps, respectively. 
N^ = B = 0. 



momentum induced by the electric field leads to an ef- 
fective magnetic field via the Dresselhaus spin-orbit cou- 
pling Eqs. ([6]) and ([7])i!- Another interesting phenomenon 
is that the spin relaxation first decreases and then in- 
creases with electric field. It is known that the electric 
field can cause two effects: hot-electron effect and center- 
of-mass drift effect i^i The hot-electron effect increases 
the inhomogeneous broadening as well as the scattering. 
The joint effects lead to enhanced (reduced) spin R&D in 
weak (strong) scattering limitj^idl when the linear Dres- 
selhaus term of the F valley is important i^ The second 
effect forces spins to precess around the induced magnetic 
field mentioned previously and hence reduces the inho- 
mogeneous broadening and the spin R&D. However, this 



To elucidate the effect of inter-valley scattering on 
spin relaxation, we compare the spin relaxation under 
four different conditions: (a) with both T-L inter-valley 
Coulomb and electron-phonon scattering, (b) with T-L 
inter-valley Coulomb scattering i?pZl only, (c) with F- 
L inter-valley electron-phonon scattering fl^r-z, o^^ly ^^'^ 
(d) without any T-L inter-valley scattering. (It is stressed 
that besides -ffp-l ^'^^ ^r-L' ^^^ ^^^ other scattering 
terms are always present.) In order to get a clear insight 
into the inter-valley scattering to the spin relaxation, we 
need to prepare enough electron density and spin polar- 
ization in the L-valleys. To do so, we first switch off the 
coherent terms at time i = (or even at i = —to, both 
yield the same result where there is no spin polarization) 
and allow the system to evolve to a spin polarized steady 
state (typically, the time tr needed to reach a steady spin 
polarization is about 6 ps after the circularly polarized 
laser pulse is applied when T — 300 K, E = 2 kV/cm, 
and B — Ni — 0). It is noted that all inter- valley scatter- 
ing should be included at this stage in order to allow F 
electrons enter into L valleys, driven by the electric field. 
Then we switch on the coherent terms to allow spin pre- 
cession and switch off the corresponding T-L inter- valley 
scattering at time t^. One should keep in mind that in 
our full calculation [such as under condition (a)], the co- 
herent terms and all the scattering terms are present at 
all the time. 

In Fig. [5] we plot time evolution of |Pa| against t — t^ 
starting from the moment when the spin relaxation be- 
gins, tr is taken as for condition (a). In the calcu- 
lation, T = 300 K, £; = 2 kV/cm and B ^ N, ^ 0. 
From Fig. [5{a), one can see that in the real situation, 
F and L valleys share almost the identical damping rate 
and precession period in spin precession, although the 
initial spin polarization of L valleys is very small. How- 
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FIG. 5: (Color online) Spin polarization \Px\ vs. t — tr under 
electric field E — 2 kV/cm. Solid curve: F valley; dashed 
curve: L valley; and chain curve: the total, (a): with both 
^r-L ^id -ffpZl; (b): without H^Z^', (c): without H^Z'l', and 
(d): without H^Zl and H^Zl- T = 300 K and Ni = B = 0. 
(e) : The initial time evolution of spin polarization of L valleys. 



ever, when H'^ZFl is absent, huge differences between spin 
relaxations in F and L valleys appear: Despite whether 
-ffp-L is present [Fig.IHb)] or absent [Fig. [Sljd)], the spin 
relaxation in L valleys is much faster than that in F val- 
ley, regardless of much larger initial spin polarization now 
gained by the L valleys due to our controlling trick pre- 
sented above [it is shown in Fig. [5l^e) that Pl (0) can even 
be higher than 5 %]. Furthermore, there is no spin pre- 
cession in L valleys anymore. It is further noted that the 
F-L inter-valley electron-phonon scattering is the domi- 
nant inter- valley scattering. This can be seen from Fig. 
[SJc) that in the presence of -ffp-L ^'^^ ^"^ ^^^ absence 
of H, 



same damping rate and precession period. This is be- 
cause when the electric field is relatively low (here E = 2 
kV/cm), only a small number of electrons distribute in 
L valleys [see the inset of Fig. [SJ a)]. Therefore the inter- 
valley Coulomb scattering is very weak. 



& 10 ■ 




FIG. 6: (Color online) Spin relaxation time r vs. electric field 



E. Solid curve: with both H^ 



_\ and H^Jj^ [case (a) 



Dashed 
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the spin polarization of L valleys shows the 



curve: without H'^J'j^ [case (b)]; Chain curve: without H^_'i^ 
[case (c)]. T = 300 K and iVi = B = 0. 

The above features can be well understood with respect 
to the different spin-orbit couplings in different valleys, in 
conjunction with the strong intcr-valley scattering. The 
spin-orbit coupling strength of L valleys is about 16 times 
as large as that of F valley4^ This strong inhomogeneous 
broadening leads to a fast spin decay in L valleys if they 
are isolated from F valley [see the dashed curve in Fig. 
[51Jd), which has a corresponding spin relaxation time in 
the order of 0.1 ps]. Moreover, the fast spin relaxation 
of L valleys also manifests itself in the quick initial spin 
decay after turning on the coherent terms in Fig. [n^b-d), 
as shown in Fig.[5]Je). It is therefore understood that due 
to the efficient inter- valley exchange of electrons caused 
by the inter-valley electron-phonon scattering as well as 
the large spin-orbit coupling in L valleys, L valleys play 
the role of a "drain" of spin polarization and enhance spin 
relaxation of the total system effectively [comparing the 
chain curves in Fig. [51^ a) and Fig. [5l^b)]. Meanwhile, the 
strong inter-valley electron-phonon scattering causes the 
spin evolutions of different valleys to be almost identical. 

To have a further insight into the effect of inter- 
valley scattering on spin relaxation under different elec- 
tric fields, the spin relaxation time r against E is shown 
in Fig. [6] under conditions (a), (b) and (c) defined above. 
In the figure, the enhancement of spin relaxation caused 
by L valleys via intcr-valley electron-phonon scattering is 
obvious by comparing the solid and dashed curves. The 
effect of inter-valley Coulomb scattering becomes visi- 
ble when E is high (> 2 kV/cm), leading to a marginal 
decrease of the spin relaxation. That is because in the 
strong scattering limit, adding new scattering tends to 



reduce the spin relaxation by its counter effect on inho- 
mogeneous broadening^ 
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FIG. 7: (Color online) (a): Temporal evolution of |Pa1 with 
electric field E = 1 kV/cm and magnetic field B = 4 T. (b): 
Spin precession frequency uip against electric field E when 
static magnetic field B = T (dashed curve), —4 T (dotted 
curve) and 4 T (chain curve). Triangular dots: spin precession 
frequency under electric field predicted by Eq.[T2]when B = 
T. Solid curve: calculated effective g-factor g* under different 
electric fields (note the scale of this curve is on the right hand 
side of the frame). 



C. Spin precession 

Now we turn to investigate the property of spin pre- 
cession. It is known that in the presence of electric field, 
spins can precess under the electric-field-induced effective 
magnetic fieldi^i. When a static magnetic field is further 
applied, spins then precess around the combined effec- 
tive magnetic field. We define the combined effective 
magnetic field term as 0(B) = JIq + 55*MsB. Here JIq 
is the induced magnetic field term, which is dependent 
on the electric field, g* is the effective g-factor. The 
corresponding spin precession frequency around ri(B) is 
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(11) 



From Eqs. ([6]) and ([7]), one finds immediately that for F 
valley, JIq is along the x-axis when the electric field is 
along the a;-axis (note that (fcp^) is larger than the mean 
value of fcp due to the strong confinement along the z- 

axis), while for L valleys, {JIq'} share the same magni- 
tude but direct along (0, ±1, ±1) directions. Therefore 
there is no net induced magnetic field under electric field 
in L- valleys, i.e., ^ JIq* = and flo of the whole sys- 
tem comes from Oq . A direct inferrer is that no spin 
precession will appear for electrons in L valleys under 
electric field only. However, as shown in Fig.[S]Ja), due to 
the strong F-L inter-valley scattering, spin precession in 
L valleys under electric field only is observed and shares 
the same precession period as the F valley. One may 
also expect that when a static magnetic field is applied, 
electron spins at L valleys may have a distinct spin pre- 
cession frequency as the g-factor of L valley is three times 
as large as that of F valley. Nevertheless, again this is 
not the case. Due to both the strong inhomogeneous 
broadening (i.e., fast spin relaxation in L valleys) and the 
intensive exchange of electrons between F and L valleys 
due to the strong inter-valley scattering, L valleys still 
present the same precession period as that of F valley, as 
shown in Fig. [TJa). This is because electron spins in L 
valleys can hardly show any spin precession characterized 
by gLfiBB/2 without being depolarized from the strong 
inhomogeneous broadening-induced interference decay as 
well as being scattered back to F valley due to the strong 
inter-valley scattering. (In fact, the spin precession time 
characterized by ql^bB /2 is 13.2 ps, the spin precession 
time under the momentum-dependent effective magnetic 
field from spin-orbit coupling in L valleys is around 0.2 
ps and the L-V inter-valley electron-phonon scattering 
time calculated using Fermi's golden rule is about 0.4 
ps when T = 300 K.) Therefore, the spin polarization 
in L valleys is just "pumped " from the F valley due to 
the strong inter- valley scattering. Consequently, the spin 
precession of the total system is determined by F valley. 
In the following, we prove this by exploring the effective 
g-factor of the total system. 

The frequency of spin precession ujp as function of elec- 
tric filed E is plotted in Fig. [7]Jb), in both the presence 
and absence of the applied magnetic field in Voigt con- 
figuration. It is shown that when the magnetic field is 
absent or is applied along — x-axis direction, u)p increases 
monotonically with the electric field. Nevertheless, when 
the magnetic field is along x-axis, ujp first decreases to 
zero and then increases again with E. The monotonic 
increase in the absence of the magnetic field is due to the 
increase of the electric-field-induced magnetic field. In 
fact, when the spin polarization is small, in low electric 
field regime the induced magnetic field term ilo can be 



roughly estimated by Eq. (|lip of F valley, i.e., 



nn 



7rnpVr 

h3 



£f 



2(l-e-^-fA-B^r) 2mra2 



(12) 



The spin precession frequency calculated using Eq. (I12|) 
in the absence of magnetic field is plotted as triangular 
dots in Fig. [7{b). It is seen that they coincide pretty 
well with the full calculation at low electric field regime. 
As the induced magnetic field is along — a;-axis, conse- 
quently an applied static magnetic field along — x-axis 
increases the spin precession frequency universally by 
a constant number (comparing the dashed and dotted 
curves in the figure). However, when the applied mag- 
netic field is along x-axis which is in opposite direc- 
tion to the induced magnetic field, the induced mag- 
netic field docs not overcome the static magnetic field 
until the electric field is larger than 2 kV/cm. By uti- 
lizing the relation -\g* ^jlbB\b=4: = ^(0) - ^(4) (or 
\g* lJLBB\B=-i = f^(-4) - rj(0)), one can get g* under 
different electric fields. Simple calculations show that g* 
changes little with electric field and is around the value 
when the electric field is abscent, i.e., g* ^ gr = —0.44 
[see the solid curve in Fig. [7](b)], suggesting that the spin 
precession properties of the total system mainly come 
from electrons of F valley. 



D. Effect of temperature, magnetic field, and 
impurity on spin relaxation 

In this subsection we discuss the effect of temperature, 
magnetic field, and impurity on spin relaxation. Spin re- 
laxation time T against electric filed E at different tem- 
peratures T = 300, 200 and 100 K with and without 
magnetic field/impurities are shown in Fig. [51 

From the figure one can see that in each case, the peak 
in T-E curve (first discussed in Sec. IIIB) is more pro- 
nounced and shows up at a relatively higher electric field 
when T is decreased. This is due to two facts: First, the 
hot-electron effect mentioned previously in Sec. Ill A is 
more important when T is low.^^ Second, when temper- 
ature is low, the F-L inter- valley exchange of electrons in 
low electric field regime is marginal due to the low elec- 
tron density in L valleys [see the inset of Fig. [51^ a)]. Thus 
the enhancement of spin relaxation by L valleys does not 
come into effect markedly until the electric field is high 
enough. 

It is also shown from the figure that in the low electric 
field regime {E ^ 2 kV/cm), the spin relaxation can be 
reduced by increasing temperature and/or magnetic field 
and/or impurity density. The reduction of the spin re- 
laxation due to higher temperature and impurity density 
corresponds to the hot-electron effect in strong scattering 
limit — the increase of temperature and/or adding more 
impurities strengthen the scattering, thus reduce the spin 
R&D. The main effect of static magnetic field in the low 
electric field regime is to reduce the spin R&D by forcing 



& 10 



B 10 




(c)B= 


4 T, Ni=0.5N 






<^--»- 




T=300 K 
200 K 




s 


-■*\ 


-■•■-- 






100 K 


-■•■■- 




\ \ 








X"* 






• 






■ 




^_ X, 








^~~~^- 










"--^ 





3 4 5 

E (kV/cm) 

30 

25 
20 

(XI 

a 15 

10 

5 



01234567 

E (kV/cm) 

FIG. 8; (Color online) Spin relaxation time r against electric 
filed E when temperature r=300 (solid curves), 200 (dashed 
curves) and 100 K (chain curves) under different conditions: 
(a) B = A^i = 0; (b) _B = 4 T and iVi=0; and (c) B = 4 T 
and Ni — 0.5N. Inset of (a): fraction of electrons in L valleys 
against electric field. 



spins to undergo a Larmor precession around it and thus 
reducing inhomogeneous broadeningjiS 

Finally one notices from the figure that when E is high 
enough, the effects of lattice temperature, impurity den- 
sity, and magnetic field are marginal. In this regime, 
the inter-valley electron-phonon scattering is consider- 
ably strengthened due to the high hot-electron temper- 
atures, and thus dominates and becomes insensitive to 
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the background temperature. Thus the changes in lat- 
tice temperature and other kinds of scattering cause httle 
discrepancy in spin relaxation. Meanwhile, the inhomo- 
geneous broadening is highly enhanced as well in this 
high electric field regime. Therefore the finite magnetic 
field has little restriction on inhomogeneous broadening 
and affects spin relaxation marginally. Consequently, the 
spin relaxation under high electric field tends to be iden- 
tical and insensitive to the lattice temperature, magnetic 
field and impurity density. 



IV. CONCLUSION 

In conclusion, we investigate the multi-valley spin re- 
laxation in n-type GaAs QWs in the presence of in-plane 
electric field under high temperature (^ 100 K), by ap- 
plying the KSBE approach. Due to the small well width 
in this investigation, electrons at T and L valleys de- 
termine the properties of spin dynamics at high electric 
field. We discuss the effect of inter- valley scattering on 
spin relaxation, the spin precession properties, and the 
spin relaxation under different conditions, such as lattice 
temperatures, magnetic fields and impurity densities. 

First, the negative differential electric conductance, 
which is essential for the charge Gunn effect, is repro- 
duced. Meanwhile, the hot-electron temperature is also 
investigated, showing that T valley has two distinct tem- 
peratures. One is for electrons in the energy regime which 
overlaps with that of L valleys, and shows the same hot- 
electron temperature of those in L valleys (Tj^). The 
other is for electrons in the lower energy regime of T 
valley (Tr). It is shown that Tr > Tl- The identical 
electron temperature in the overlapped energy regime of 
r and L valleys is due to the strong inter-valley scatter- 
ing whereas Tr > Tl results from the smaller effective 
mass in F valley. 

The spin relaxation time t first increases and then de- 
creases with the electric field, especially at low tempera- 
tures. The first increase of r in low electric field regime 
is due to the hot-electron effect in the strong scattering 
limit. Whereas the following decrease of r in the high 
electric field regime is mainly due to the fast spin relax- 
ation in L valleys (originating from the much stronger 
spin-orbit coupling at L valleys) and strong T-L inter- 
valley scattering. In fact, L valleys are shown to play the 
role of a "drain" of the total spin polarization. It is also 
shown that the T-L inter-valley electron-phonon scatter- 
ing is responsible for the inter- valley scattering. The con- 
tribution of T-L inter-valley electron-electron Coulomb 
scattering is marginal as it does not cause any electron 
exchange between different valleys. Moreover, due to the 
strong F-L inter-valley scattering as well as the strong 
inhomogeneous broadening in L valleys, electrons in L 
valleys possess almost the same spin relaxation rate and 
spin precession period as those in F valley. It is further 
shown that this feature still holds when a static magnetic 
field is applied. Actually, in our investigation {B < 4 
T), the spin precession frequency is determined by elec- 



trons at F valley, despite the fact that gL is much larger 
in magnitude than gr- This is because the net electric- 
field-induced magnetic field from L valleys is zero, and 
the effect of static magnetic field in L valleys can not 
surpass the much stronger inhomogeneous broadening as 
well as frequent exchange of electrons from F valley. 

We also investigate the effects of temperature, mag- 
netic field and impurity on spin relaxation. When the 
electric field is not too high, the spin relaxation can be 
reduced by increasing temperature, magnetic field, and 
impurity density. However, when the electric field is high 
enough, spin relaxation time is dominated by the F-L 
inter-valley electron-phonon scattering and becomes in- 
sensitive to temperature, impurity density and magnetic 
field. 

Finally we remark on the possibility of the "spin Gunn 
effect" in GaAs QWs. First, as pointed out in our 
previous work 4^ the feasibility of drift-diffusion model 
used in Ref. [SJ is valid only in system without Dressel- 
haus/Rashba spin-orbit coupling. For Zinc-blende semi- 
conductors, even without magnetic field, electron spins 
can flip due to the Drcsselhaus/Rashba spin-orbit cou- 
pling, let alone an effective magnetic field induced by 
electric field via the same spin-orbit coupling. There- 
fore, although with the spin-dependent mobilities, elec- 
trons with opposite spins can experience frequent spin- 
flip. Second, our calculations show that in high electric 
fleld regime, spin relaxation time decreases effectively 
with electric field. This is due to the hot-electron ef- 
fect, together with the fast spin relaxation in the upper 
valleys. Therefore, the condition for "spin Gunn effect" 
proposed in Ref. |34| can hardly be realized due to the 
short spin relaxation time in GaAs QWs. 
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APPENDIX A: THE COHERENT AND 
SCATTERING TERMS OF THE KSBES 

The coherent terms are 

PkJcoh = -i[^\{'kx)-fT-^Vxx,qPk^-q,Pk^] , (Al) 

q 

with riA(kA) = hA(kA) + iffAA^sB and V^A.q represent- 
ing the intra- valley Coulomb scattering matrix element. 
[A, B] stands for the commutator of A and B. The scat- 
tering terms are 



Pk;, 



Uat.xy = -{S^^A>^<) - StU<^>) + S^^,,i>,<y ^ S^^,,i<,>)n , 



(A2) 



where 



q 

+ ^Y1 ^AA',qPk;,-qPk;,Tr[p>, p< _q](5(ek;,-q - £k;, + £k;, - £k'^,-q) 



q.k'w 



'^E^AA',qPk;_,PkJ^Al.''5(£k. -£k'^, -f^AA') + ^A''A''5(£k. -£k'^, +f^AA')] 



(A3) 



r 



Here p^ = pk and p^ = 1 - pk- A^aa' = ^a'a = 
[exp{fiCt\\> /ksT) — 1]~^ is the Boson distribution 
of phonons with frequency f^AA'- N^y = Nxx' and 
1 + N\\i. Mxx'.q is the matrix element of 



^A>A 



K 






electron-phonon scattering with q = k^ — k'^, 
being the two-dimensional phonon wave vector. The 
intra-valley electron-acoustic-phonon scattering is ne- 
glected in the present investigation due to the high 
temperature. The matrix elements of the intra-valley 
electron-longitudinal-optical (LO) phonon scattering 
in r valley and the intra-valley optical-phonon de- 
formation potential scattering in L valleys are given 



by M^r q 
Mf r „ = 






V^OQ 






K 







Irr{iqz)\'^ and 
^, respectively!^ For 



the inter-valley electron-phonon scattering, selection 
rules applied to a cubic zinc-blende structure of III-V 
semiconductors show that the LO and longitudinal- 
acoustic (LA) phonons can assist the inter-valley 
transitions;^ However, it was further shown later 
that the selection rules do not apply in the high energy 
regime where the inter- valley transfer can happen j^^'^^i^^ 
Therefore, all phonon branches including transverse op- 
tical (TO) and transverse acoustic (TA) phonons can 
contribute to the inter-valley scatteringi^ Nevertheless, 
the electron- TA phonon scattering is very weak and can 
be neglectedf^iS The LO, TO and LA phonons have 
comparable energies?^ thus the inter-valley scattering 
can be considered by grouping together the scatter- 
ing assisted by these phonons and using the average 
phonon energy and the total coupling constant.— Here 
the two-dimensional phonon wave vector q is mainly 
determined by the large component k° — k°, , thus 
the average phonon energy and the total coupling 
constant are approximated to be fixed. Therefore, 



we have M^j^, 






j:. 



noi 



Ti 



i^<lz)\ 



for T-L inter-valley electron-phonon scattering, and 

KDl ^ 

^n^L^.o, = Eg, idQ.l^i,V uLAnz)Y for L-L inter-valley 
electron-phonon scattering. Here d = 5.36 g/cm'^ is the 
mass density of the crystal; kq = 12.9 and Kqo = 10.8 
are the relative static and high-frequency dielectric 



constants respectively; cq is the vacuum dielectric 
constant. The phonon energies are MIyt ~ 35.4 meV, 
hVlL,u = 34.3 meV,.^ hVlTu = ^^L.r = 20.8 meV,-^ 
and HUl-l- = 29.0 meV^ The deformation poten- 
tials are Dl,l. =0.3x10^ eV/cm,^^ Dtl, = 1-1 x 10^ 
eV/cm^i^ 'and Dl,Lj = 1-0 x lo'^ eV/cm^iS 

VxA'.q = 'Eq, ^^^„(qi\qi + ^-l) h\{iqz)Il,y{iqz) 

Eq, eo.o{q-iqi+^-) \^^^^^1^-)\^ ^^ ^^"^ intra-vallcy 
(inter-valley) Coulomb scattering matrix el- 
ement when A = A' (A ^ A')- C/|a q = 
EqAZieVieoKoiq^ + ql + n')]}^\hx{^qz)\'' ^ is the 
intra-valley electron-impurity scattering potential, with 
Zi — 1 being the charge number of the impurity. 
Ni is the impurity density, k? — Ne'^ /{asoHokBT) 
is the screening constant. The form factor 



with y = a{qz 



Kl. 



Kl, 



)/2. 



y2(^y2 



APPENDIX B: NUMERICAL SCHEME 

We extend the numerical scheme for solving the KSBEs 
at F valley only, shown in detail in Ref. [ll|, to the case 
with both F and L valleys. To do so, five circular zones 
centered at the bottom of each valley in momentum space 
are set up. We divide the truncated two-dimensional 
momentum space in F valley into Nr x M control regions, 
and that in Li {i—1, ..., 4) valley into N]^ x M control 
regions. The relative k^-grid point in A valley is chosen 
to be the center of the region. 



kn.m 



n 



(cos6'„,sin6'™). 



(Bl) 



with E^ ^ {n + 0.5)AE and e„, = mM (n=0, 1,- • • , 
iV^t - 1 and m=G, 1, • • • , M - 1). Here 7Vf„t = Nt 
for F valley and N^^^ = Nl for L valleys. The trun- 
cation energy is E^^ _ . Due to the effective T-L en- 



^w.-i 



ergy difference, E^ -^ 

Nl are chosen to satisfy Ny — Nl 

integer part of ^^). The energy partition AE is set 



ErL- Thus N-p and 
f^i^l ([fffcl is the 



lAEi 



lAEi 
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as Mlrr/nLO with tilo being an integer, and the angle 
partition is A0 — 2ti/M . The coherent terms and the 
scattering terms of electron-impurity scattering can be 
discreted directly. By virtue of the chosen energy par- 
tition, inravalley electron-phonon scattering in T valley 
can be discreted immediately as well. For intra-valley 
electron-phonon scattering in L valleys and inter-valley 



electron-phonon scattering, one expects to have the en- 
ergy difference be the integer multiplication of A£'. This 
is difficult to satisfy in general. However, with an opti- 
mal value of TT-LOj one can approximately achieve above 
requirement. The remainin g n umerical scheme is all the 
same as those given in Ref. [ill 
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